The electro-magnetic pulse produced by a nuclear bomb 

explosion high above the atmosphere. J. Steinberger Sept. 2006 



Introduction. 



This article discusses the electromagnetic pulse produced by a high, exoatmospheric nuclear 
explosion, commonly abbreviated HEMP. The gamma rays produced in the nuclear explosion above 
the atmosphere produce electrons in the atmosphere by Compton scattering. The electron current, is 
initially downward. This downward current does not produce a downward propagating electro- 
magnetic field, but its deflection in the earth magnetic field produces a transverse component of the 
current, which does result in an electromagnetic wave propagating downward. This was realized in 
the early 60's, in connection with nuclear bomb tests well above the atmosphere, in particular the 
"Starfish" explosion in 1962, 1 .4 Megatons, 400 km above Johnston Island in the Pacific Ocean. The 
atmospheric test ban of 1963 has prevented further explosions of this sort. The U.S. Department of 
Defence, in unclassified reports 1) , notes that in connection with the Starfish test 30 street light circuits 
went out on the Hawaiian island Oahu, on which Honolulu is located, 1200 km from ground zero. 
The electric fields of the pulse is canonically stated 1 ' to be -50,000 volts/meter. It is often asserted 2 ' 
that such an explosion, although it would not harm the population, would disable the electric power 
and communications networks of a whole continent, and so make appropriate military response 
impossible. Given the military and political implications of this pulse, as well as the fact that no 
nuclear weapons expertise is required to calculate it, I thought it worth while to try to check these 
assertions. In the following, cgs units are used. 

HEMP produced by a high altitude nuclear explosion. 

Field equations . 

In adequate approximation, the problem of the electromagnetic wave produced below it by an 
explosion at several hundred km altitude, the surface of the earth can be taken as flat. The solution 
to this problem has been given by Karzan and Latter 3 ' and by Longmire 4 '. Here we follow the same 
physics. The relevant Maxwell's equations are 

1 / c8B y I dt = -8E X /dz; 1 / cdB x ldt = 8E y l dz. ; 
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\lcdE x ldt + E x + — j x =-8BI8z, \/cdE y /dt + E+ — j =+dB x /dz, units are cgs. 
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t and z are time and distance from the explosion, a is the conductivity due to the ionization 
produced by the Compton electrons. 
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It is useful to introduce the retarded time 3 ' 4 ' x , x=t-z/c. - — >- 



dt dr dz dz cdr 
The field equations become: 
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a} — = -— ^ + — ^; b) —^ +—(cjE x + j x ) = -—^ + — , and 
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dB x dE v 8E y 8E y 4/r . v . . dB x dB x 
a) — - = — 1 -; b) — - + — (oE +j ) = — ; 
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Adding a) and b), and subtracting a') from b'), also noting that for the plane wave moving along z, 
B } = E x , and B x = -E y , there is finally the very simple result: 
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(1) 8E i /8z + E ( + — ; ( =0 ; where i = xory. The electric field E, the current j and the 
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conductivity o are functions of z and the pulse timer . 



Gamma flux. 
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Fig. 1 . Geometry 

The assumed geometry is shown in Fig. 1 . The earth surface is taken as a plane, the 
curvature is ignored. The angle of observation with respect to the vertical is y The coordinate 
system is taken with the z direction from the ground to the burst, along the direction defined by the 
angles , and the x direction in the plane defined by z and the magnetic field. The angle of the 
magnetic field with respect to x is 0. The height of the explosion is h. In the following, for numerical 
examples, h is taken to be 500 km and to be 45°. 

The average gamma energy is ~ 1 .5 MeV. The number of prompt gamma rays produced by a 
nuclear explosion is between .3% and 1% of the energy released in the explosion 4) Here the 

upper value is taken, which corresponds to lxl0 26 gamma rays per Mton yield, with average gamma 
ray energy 1 .5 MeV, and with burst duration about 20 nano-seconds. The gamma rays Compton 
scatter in the atmosphere, transferring their energy to electrons, with a cross-section 6 ', 

cr =0.17x10 24 cm. The gamma flux as function of altitude, f (z) , is: 
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where z is taken from the ground towards the explosion, f y is the incident gamma flux, which, for a 
1Mton explosion is then: 
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-cosy 2 =3.2xl0 9 cos7 2 , for h = 500 km, 



p(z) is the atmospheric density, p(z) = .0013e^ cosr/6 Mm g I cm 3 , 
gamma flux as function of atmospheric height is then: 



and A is the Avogadro number. The 



f r (z) = 8x10 lh xcos ^exp(-45/cos7xexp-zcos//6.8fcm) . 
In the following, for convenience, zcosy /6Mm will be replaced by the dimensionless z. Then, for h = 
500km, 

(2) / =3.2xl0 9 cos 2 7exp(-45/cos/xexp-z) 
Transverse Compton electron current . 



As already noted, the down moving Compton current does not produce a down moving 
electro-magnetic wave. However, the current spirals around the earth magnetic field, thus producing 
transverse components. For a typical magnetic field of the order of .5 gauss, at 45° to the horizontal 



plane, and an average electron momentum of 1 .1 MeV/c, the Larmor radius is 5000 cm. At the 
relevant atmospheric density of ~ 2x10 ~ 5 gl cm 3 , typical electron range is of the order of 20,000 cm, ~ 
4 times the Larmor radius. The component of the current along the field is undeflected, the 
component perpendicular to the field spirals with the Larmor radius. In the following the Compton 
electron total energy (including rest energy) spectrum is approximated by a uniform (top hat) 
distribution from .5 MeV to 2 MeV, with amplitude so that the electron kinetic energy flux equals the 
gamma energy flux. The consequent transverse currents are: 

~ 2MeV ,s 

2 r (1- cos cot) 



j x (r) = ecf r (z) I X * cos sin cos 6 * — * 
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j y (t } = ecf y (z) I X * cos cos 6 * — * J dEfi 
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x is the retarded time. t(0 = (l-/?(sin # + cos 0(smm)/m))t . 
t is the time after the scattering of the gamma ray on the electron. 
A is the mean free path of the gamma rays, 

(4) ^ = ^^> 

f3 is the electron velocity divided by the light velocity c, 

£ is the average production angle of the Compton electrons, £ = 30°, 

co is the Larmour frequency, m = — P cB (s auss ^ — 

10 4 /3* pc(MeV) 

E min is the minimum energy for which the Compton electron is still in motion at pulse (retarded) time 
t . The lifetime, At, of the electron is 

^ electronrange (E - .5MeV) * .5 gem 2 MeV 1 

~ Jc ~ Pcp(z) ' 

E is the total energy of the electron, including its rest energy, and consequently, 

E mn =(7.2*10 7 ?(r)(sec)*e^ + .5) MeV 

Conductivity. 

The conductivity is due to the ionization electrons produced by the Compton electrons, ~ 
30,000 ionization electrons per MeV Compton energy, with average kinetic energy of the ionization 
electrons E ionfi = 15eV . The consequent conductivity is 

n eAv 

a = — — , where n e is the density of ionization electrons and Av x is their average velocity 

E , ' 

in the x,y direction, produced by the Field E x , y , between collisions. 

n (z) = *\.5MeV * 30, OOOionizationelectrons I MeV = 45,000^^ ; 
X X 

Av *,y =^*- E *,v* At co„. J At coii. = 7 ! An <* , ^ ; 

2 m v ion *G coU *AI\5*p(z) 

the factor V2 accounts for averaging over the time between collisions, 
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v ion = velocity of the ionization electrons; initially this is: v ion = J — * c = .0077c . 
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and o coll =6*10 16 cm 2 is the cross-section for the scattering of 15 eV electrons on nitrogen or 
oxygen atoms 7 '. The conductivity is then: 



2 

(5) a = f (z)l A * 45,000 * — * At coll ; with A? co „ the time between collisions. 

2m 

The conductivity is inversely proportional to the electron velocity v ion . This decreases with 

time, since in the collisions with the gas molecules, the electron looses, per collision, on the average, 
a fraction of its energy equal to the electron mass divided by the molecule mass, approximately 
1 .8*1 0" 5 . This is a non-negligible effect. 

1 dE. , A 

TT = -l-8*10- 5 *v ( , nCTco „.* — *p(z), 

£,„„ dt 15 

UlO'sec 'e 2 = 1/A 
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where E ion and E ion fi are the energy of the ionization electron as function of time, and its initial 
energy, that is 1 5 eV. 

Solution: E ion = E ion >0 /(1 + /7A) 2 , where A = 1.5 * 10~ 8 sec* e z , 



V- = V- n 
/on. ;o«,0 - 
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The average velocity is then: < v. >= — \dtv ion dt , where 



At Q 

15MeV * 5s * cm~ 2 MeV 1 _ s 
Ar = - — ; = 1.05*10 V is the average time of contribution of electrons to the 

.92 * c*. 0013 g* cm 3 e z 

current. The near equality of A and At seems an accident of nature. Then v ion = J6v ion Q 

A correction of the opposite sign is due to the electric field, which, between collisions, 
increases the average energy of the ionization electrons during the pulse. The average energy gain, 

per collision, is AE ion = = (eE x At coll f 12m, and dE inn /dt = AE ion /At coll , 
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y i0 n(t)=y i0 n.o^+t/t , — = — 

2/3((l + At/t ) 3 ' 2 -1) „ W/1 ^ A , 

< v ta >= v,,„, — ^ '- = v taf (1 + 1/4* A/// ) . 
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Depending on the angle 7 of the propagation being considered, as well as on the direction of 
the magnetic field, this correction increases the velocity of the ionization electrons by from 2% to as 
much as 20%. In the following, an average increase of < v ion > of 10%, from .0077*. 76 *c=.0058*c to 
.0064*c has been used. 



Finally, the conductivity used in these calculations is: 

friz) e 2 45,000 



(6) a = ^ — * * 

A, 2m .0064*c*a eo „ *A/15*p(z) 

Back reaction of the electric field on the current which produces it . 



The transverse momenta of the compton electrons constituting the transverse current are 

dp , v 

reduced by the electric and magnetic fields which they produce. — — = eE x (1 - — ) . Using the 



dt 



rough approximation for the electron velocity, Av x >y = -c 



4p, 



- * (1 - v z I c) , where c^jp 2 +(mc) 2 



lp e +(mc) 2 

is the average total energy of the Compton electrons composing the current, ~ 1 .2 MeV. The 



average reduction in the x,y component of the electron velocity is: 



<Av,, y > c eE x Jl-v z /c)A't 
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where A't = 



E y (MeV)*.5gcm- 2 MeV 1 Alg I cm : 
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and vjc = cos<§ */?(cos 2 Ocosurt + sin 2 0) . 



The effective current is: j wff = j x<y + j xMreaction ; with 
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EMP field strength and pulse duration . 
The field equations become: 



dE x 6.8km ^ 2n<j 
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With the above expressions, (6) for <r , for j x y , (7) for j x 
equations become: 

dE r 



y ,backreaction ' 



(2)for^(z), (4) for A, the field 
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with a=29.4 per unit z, b=16.4 per unit z and c=435 -J ergcm /cm 2 per unit z, assuming a magnetic field 
strength of .5 gauss, a height of 500 km, and a bomb strength of 1 Megaton 

Figure 2 shows the pulse shapes for E y and B x and figure 3 those for E x and B y , for different 
angles with respect to the vertical and consequent differences in distance from ground zero, taking 
the magnetic field angle as 45° and the gamma ray pulse duration to be 20 nanoseconds. For a 
pulse duration of 10 nsec the pulse heights are larger by 10% and the pulse length shorter by 20%, 
leaving the energy flux unchanged. 
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Figure 3. E x and B y as function of time for various observer directions y , for a 1 Mton bomb at 500 
km height and a magnetic field of .5 gauss at 45° to direction of propagation. 



Figure 4 shows the variation of pulse heights and pulse width for directly overhead explosions 
(cos y =1) at 500 km altitude, as function of the bomb strength, over a range of two orders of 
magnitude, .1 Mton to 10 Mton. It illustrates the saturation effect: for explosions above a certain 
strength, the pulse heights vary by less than a factor of two. The Pulse height is limited by the 
reactive effect of the ionization produced by the current. This limits the HEMP pulse heights which 
can be achieved to those calculated here. In particular, specialized HEMP weapons of substantially 
higher HEMP yields, as are sometimes imagined, are not possible. 

The summary table below gives the total energy flux, the pulse heights and the pulse width in 
the EMP produced by a 1 Mton explosion at 500 km, as function of the observer angle. The 
magnetic field of .5 gauss is taken to be at 45° to the direction of propagation. 



cosine of observer direction 
energy flux, ergs/cm 2 
pulse height, Kvolts/m 
pulse halfwidth, usee, 



Summary Table 
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47 
27 
.032 
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53 
20 
.065 
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Figure 4. The variation of pulse heights (left) and pulse width (right) with bomb strength (above 
ground zero, height 500 km, theta =450, geomagnetic field .5 gauss). 



Damaging effects of HEMP on electronic instruments, telecommunication and electric power 
networks 



Observation of the Starfish test on Hawaii. 

Although such a high altitude nuclear explosion would not directly endanger human life, it is 
common, in military and governmental circles to credit the HEMP with devastating effects on the 
communications and electric power networks of a whole continent, which could disable the society 
and the ability to respond to such an attack, unless these devices are especially shielded or 
protected. The effects to be expected from such electromagnetic pulses are however far from 
obvious, and difficult to calculate or to evaluate experimentally. 



July 8 1962, 1 1 PM the "Starfish" test, a 1 .4 Megaton nuclear explosion at an altitude 400 km 
above Johnston island, could be "observed" from the Hawaiian islands, 1200 km to the west. This 
event had been announced previously and according to the newspapers, the population joyfully 
anticipated it. According to the Honolulu Star-Bulletin of the next two days 8) , the sky lit up bright 
blue, as in daylight, for some seconds, streetlights on Ferdinand St. In Manoa and Kawainui St. In 
Kailua went out the instant the bomb went off, and "police were sent off to a South St. warehouse, 
when a burglar alarm started ringing at the time of the blast".. Fuses blew out in Kasmuki, Kahala, 
Kaluji, Maili, Wajanae, Makaha, Wahsaura, Kailua and Sandy Beach, and the Federal Aviation 
Agency reported occasional radio communication blackouts over a period of two hours. 

These reports, of Honolulu street lights going out, burglar alarms sounding, and circuit 
breakers popping are frequently cited as evidence for the danger posed by the HEMP, but such 
things happen also without nuclear bomb explosions, and the connection cited with the nuclear 
explosion of these incidences are far from clear. The outstanding fact seems to be that life went on 
normally in Hawaii, during and after the Starfish test, the Honolulu Star-Bulletin appeared normally 
the next day. In response to such an article citing these alleged Starfish effects 9) , there are answers 
by Alan S. Lloyd, Manager, Consumer Services, Hawaiian Electric Company 10) stating categorically 
that "there was absolutely no adverse effect on our utility power and communications system on the 
island of Oahu where the city of Honolulu is located", and by E.C. Schoen, Chief Engineer, Hawaiian 
Telephone Company, that, in connection with these nuclear bomb tests "No one, that I have talked 
to, can recall any increase in troubles associated with telephone sets or any other type of difficulty" 

Damage to microelectronics. 

One of the often cited dangers of electro-magnetic pulses is damage to microelectronics, a 
very important element of our modern society. Given the modest sizes of electronic instruments, 
such as computers radios, telephones, the effects of pulses such as those calculated above are 
relatively easy to test in the laboratory, and some colleagues here at CERN, the European centre for 
particle physics, in their spare time, are trying to study the effects of electric and magnetic pulses with 
waveforms as expected in the HEMP, on several highly sophisticated microelectronic circuits as well 
as some cell phones. Until now they have not succeeded in damaging a circuit, with voltages up to 
130 kV/m, although at the higher voltages the processes in progress could be disturbed 1 . This 
work is continuing. 

Damage to telecommunication and electric power networks . 

As noted above, the Hawaiian systems were not disrupted by the Starfish (or other) high 
altitude nuclear test. I am not capable to estimate the damaging effects to be expected. It should be 
noted that although it is sometimes asserted that the voltages could propagate and build up on the 
network lines for hundreds of kilometres, and therefore become extremely large, it is my 
understanding that this is false, that the electric field cannot build up over greater distances than the 
order of the coherence length, that is, the pulse length, 50 nanosec*c = 15m. 

Strong electromagnetic pulses are produced by lightning currents. The voltages are typically 
greater than the EMP voltages calculated above, for distances within a few hundred meters of the 
current, and since they are ~ 100 microseconds long, compared to the 50 nanoseconds of the 
HEMP, their energy flux is correspondingly greater. Lightning strokes within several hundred meter 
distances occur typically once per several years, yet I have not succeeded in finding clear evidence 
of damage caused by lightning electromagnetic fields, as distinguished from the effects produced by 
the lightning currents, to telecommunication networks, power distribution networks, or 
microelectronics, 



Conclusion on damages . 



It is not clear that spending billions, as some propose 13) to congress, on hardening of 
telecommunication, electric power networks, and public microelectronics to "protect" against a 
possible HEMP, is in the public interest 
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